Pressure-induced hole doping in the Hg based cuprates 
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We investigate the electronic structure and the hole content in the copper-oxygen planes of Hg 
based high T c cuprates for one to four CuC>2 layers and hydrostatic pressures up to 15 GPa. We 
find that with the pressure- induced additional number of holes of the order of 0.05e the density of 
states at the Fermi level changes approximately by a factor of 2. At the same time the saddle point 
is moved to the Fermi level accompanied by an enhanced k z dispersion. This finding explains the 
pressure behavior of T c and leads to the conclusion that the applicability of the van Hove scenario 
is restricted. By comparison with experiment, we estimate the coupling constant to be of the order 
of 1, ruling out the weak coupling limit. 

PACS numbers: 74.25.Jb, 74.72. J, 74.62. Fj, 71.15.Mb 



In the high temperature superconductors, the critical 
temperature strongly depends on the number n of CuC>2 
layers per unit cell. In particular, an increase of the 
transition temperature with increasing n has been ver- 
ified within the different cuprate families like the Bi, Tl, 
or Hg based compounds. Although it seemed tempting 
to add more and more such CuC>2 layers to the crystal 
structure, there is a limit in that for n > 3 the super- 
conducting properties become worse, i.e. T c decreases 
again. Besides the composition, pressure is an impor- 
tant tool to enhance the transition temperature. This 
effect is most pronounced for the Hg based cuprates 0, 
which still hold the world record in T c being 97 K for 
n=l and ambient pressure, and even 164 K for n=3 at 
30 GPa 0, El ■ Therefore, this system is best suited to 
simultaneously study both effects. 

In this Letter, we report on the electronic properties 
of HgBa2Ca,i_iCu Tl 02n+2 as a function of composition, 
i.e. n = 1,2,3,4, and hydrostatic pressure, investigated 
by first-principles calculations. Special emphasis is lead 
to the question how the charge carriers are redistributed 
leading to an increase of the hole concentration in the 
CuC>2 layers, which are responsible for the superconduct- 
ing current. At the same time the density of states (DOS) 
at the Fermi level is considered to shed light on the role 
of the van Hove sing ularity (vHS) for the supercon- 
ductivity and to estimate the coupling constant. 

All calculations, which are based on density func- 
tional theory (DFT), have been carried out with the 
full-potential linearized augmented planewave (LAPW) 
method utilizing the WIEN2k code The atomic- like 
basis set used in the atomic spheres around the nuclear 
positions [fl Q allows for an analysis of the electronic 
charge in terms of atomic origin, and for a decomposi- 
tion of the charge according to its orbital momentum I. 
Therefore, the decrease of the corresponding occupation 
number can be interpreted as the creation of a certain 



amount of holes of the particular character. Below, the 
Cu02 layers will be analyzed that way. The most im- 
portant orbitals are Cu-d K 2_, y 2 and 0-p x forming pd — a 
bonds. Note that the 0-p y orbital of the second oxygen 
neighbor is equivalent to 0-p x orbital of the other oxy- 
gen due to symmetry. For all our calculations, the crystal 
structures have been relaxed in terms of lattice parame- 
ters and atomic positions as determined and discussed in 
detail in Ref. @ This procedure makes the calculations 
not to rely on experimental crystalline data. The latter is 
important for studying the pure dependence of the elec- 
tronic properties on pressure or n, since the measured 
samples usually include some excess oxygen and can be 
strongly Hg deficient Figure Q depicts the number of 
holes created in the Cu-d x 2 _ y i orbitals when pressure is 
applied. Before discussing the details, it should be noted 
that for all compounds and pressures the atomic sphere 
radii, chosen as 2.1, 2.3, 1.9, 1.2, and 2.0 a.u. for Hg, 
Ba, Cu, O, and Ca, respectively, in which the electronic 
charge is analyzed, has been kept constant. This means 
that a bigger fraction of the crystal volume and hence 
electronic charge is lying inside the sphere when pressure 
is applied / increased. Therefore, the number of pressure- 
induced holes can be regarded to be somewhat underesti- 
mated by our pressure calculations. Nevertheless, it can 
be stated that there is a pronounced pressure effect on 
the orbital occupation numbers of copper, whereas there 
is hardly any changes in the corresponding oxygen 0-p x 
charges 8j. This is in contrast to the doping-induced 
holes p.0] which are found not only at the copper sites, 
but also, to some extent on the oxygen positions as well. 

Let us first focus on the partial charges at ambi- 
ent pressure. When changing the composition, the Cu- 
d x 2_ y 2 occupation number changes from 1.4336 for n = 1 
to 1.4414 for n = 2, 1.4421 (1.4386) for n = 3, and 
1.4413 (1.4380) for n = 4. The numbers in brackets re- 
fer to the inequivalent copper positions occurring for the 
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FIG. 1: Number of pressure-induced holes (in units of elec- 
trons) for the Hg based cuprates, Hgl201, Hgl212, Hgl223, 
and Hgl234 for pressures up to 15 GPa. The values are taken 
with respect to ambient pressure for each material. The filled 
circles correspond to the total number of holes created per 
unit cell, whereas the open symbols (triangles for innermost 
and squares for outer layers) reflect the situation in a partic- 
ular layer. The symbols indicate the calculated data points, 
while the lines serve as guide to the eye. 



three and four layer compounds and indicate an inho- 
mogeneous hole distribution that concerns the different 
CuC>2 layers within one unit cell. Upon applying pres- 
sure, the number of holes increases with a parabolic-like 
shape which at 15 GPa is close to its maximum for n up to 
three. Only for n—A, there is already a slight decrease in 
the highest pressure range. An important fact is that the 
amount of pressure- induced charges per site decreases for 
higher n. In particular, the innermost layers (marked by 
triangles in Fig. do not gain more holes with pressure 
beyond n = 3. Only the total amount of holes per cell 
created for a certain pressure value is similar in all the 
compounds, i.e. when the contributions from different 
layers are summed up weighted by their multiplicity. 

All these findings are fully concomitant with the super- 
conducting transition temperature [(|: Experimentally, 
the enhancement of T c for n up to 3 with pressure is 
approximately 15-20 K between and 15 GPa, i.e. it 
is nearly independent of n. Only for n—A the curve is 
somewhat flatter. Exactly the same trend is found for 
the number of pressure-induced holes. The inhomoge- 
neous charge distribution among the different layers also 
provides an explanation for the limited T c since the inner 
layers do not gain any more holes by further insertion of 



Cu02 planes. 

The fact that T c scales with the number of charge carri- 
ers does, however, not provide an explanation for the ori- 
gin of superconductivity. Although the amount of holes 
in the cuprate planes is an important quantity, its knowl- 
edge does neither lead to a conclusion about the pairing 
mechanism, nor about the coupling strength to the medi- 
ating quasiparticles. Therefore, the question arises what 
the role of the bands and the density of states around the 
Fermi level is. In this context, the van Hove scenario has 
been frequently invoked. Indeed, the most precise calcu- 
lations for the bandstructure of YBa2Cu307 0] exhibit 
a saddle point (usually referred to as a vHS) pinned to 
the Fermi level, and the same is true for La2_yBa T Cu04 
[l2| and HgBa2Cu04+5 at optimal doping |1Q| . 

Here, we will assess this possibility for the Hg based 
compounds under pressure by performing selfconsistent 
calculations for all pressure values to obtain the corre- 
sponding electronic structures. Thereby we observe two 
effects: First, upon applying pressure, the dispersion in 
the k z direction is strongly increased. The band width of 
the pd — a band along XR is changed from 0.12 to 0.29 
eV between and 10 GPa which corresponds to a factor 
of 2.5. Second, at the same time the ratio of hopping ma- 
trix elements t' jt decreases by 6%, where t 1 and t denote 
the in-plane matrix elements of 0-0 and Cu-0 hopping 
when mapping the band structure to a simple tight bind- 
ing model. These two effects lead to an asymmetry in 
the location of the vHS with respect to the Fermi level 
for different k z values. This can be seen in Fig. |2 where 
the pressure dependence of the Hgl201 band structure 
around Ep along the two high symmetry lines TX and 
ZR is depicted. Upon applying pressure, the shoulder 
between Z and R clearly moves up toward the Fermi 
level intersecting it at around 10 GPa. Along TX the 
band is pushed up as well thereby getting steeper, but 
staying below Ep in the entire pressure range. As a re- 
sult, the vHS is broadened when it arrives at the Fermi 
level. Thus, we have found the vHS to be pinned to the 
Fermi level within a certain region of the Brillouin zone 
(between TX and ZR) at high pressures which could be 
related to the high and rather weakly pressure-dependent 
T c in a pressure range between 10 and 20 GPa. 

These findings also show up in the density of states de- 
picted in Fig. For Hgl201, the shoulder in the DOS is 
located at -0.5 eV but moving up toward Ep with pres- 
sure, where the density of states at the Fermi level Ao is 
only slowly varying for higher pressure values (see Fig. 
0J. For n=3, the situation is the same, while for n=2 
and 4 the vHS crosses Ep between ZR and stays below 
Ep along TX. Concerning the magnitude of the DOS 
the following is interesting to note: Ng is getting larger 
with the number of layers, where the pressure-induced 
increase is nearly the same for n=l, 2, and 3. The rea- 
son for the latter fact is that with additional Cu02 lay- 
ers the corresponding bands are not fully degenerate but 
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FIG. 2: Band structure of Hgl201 along the high symmetry 
lines FX and ZR for different pressure values. The black 
filled circles mark the calculated points, where the dotted lines 
serve as guide to the eye. The Fermi level is set to zero and 
indicated by a horizontal line, and the energy is given in eV. 

exhibit a certain splitting which does not allow all the 
corresponding shoulders (vHSs) to sit at Ep at the same 
time. Thus, the pressure-effect on the DOS is indepen- 
dent of the number of layers, at least up to n=3. This 
trend goes hand in hand with the number of holes as well 
as with the critical temperature. Only for rt=4, the pres- 
sure effect is less pronounced, since at ambient pressure, 
Nq does not fall into a local minimum as is the case for 
the other compounds. 

We want to point out that the selfconsistent treat- 
ment is crucial for understanding the changes of the 
band structure under pressure. For example, when using 
the rigid band approximation to estimate the pressure- 
induced number of holes necessary to push the vHS to 
the Fermi level, as it was proposed first by Novikov et al. 
[l3|. this quantity being of the order of 0.5e is overesti- 
mated by about an order of magnitude (see Fig[5J). 

Knowing on the one hand, the pressure dependence 
of N for the four compounds as a result of our calcu- 
lations, and the corresponding experimentally measured 
T c s on the other hand, one can try to estimate the effec- 
tive dimensionless coupling constant A = VNq, assum- 
ing a BCS-like behavior valid for the weak coupling for 
any type of the intermediate boson, T c ~ lob exp(— 1/A). 



FIG. 3: Densities of states TV in states per eV and unit cell 
of Hgl201, Hgl212, Hgl223, and Hgl234 for pressure values 
of 0, 5, 10 and 15 GPa. 



Here, lob is the frequency of the mediating boson, and V 
characterizes the coupling strength. It should be noted 
here that the DOS does not necessarily have to be as- 
sumed to be constant near the Fermi level on the cob 
scale, but it turned out to be approximately the case (Fig. 
|3J and thus can be used for such an estimation. There 
are, however, several problems which hamper a quantita- 
tive comparison between theory and experiment, (i) The 
main concern comes from the fact that T c has been mea- 
sured on samples which include a doping concentration 
close to optimal doping. This doping could be considered 
as an unknown chemical pressure. On the contrary, the 
calculations have been performed for the undoped mate- 
rial, (ii) There are big uncertainties in the experimental 
data regarding defects and non-stoichiometry of the sam- 
ples, (iii) When dealing with compounds with more than 
one layers, it is unclear how their contributions enter T c . 

Despite all these uncertainties, conclusions are never- 
theless possible. First, from all calculations performed to 
investigate doping as mentioned before for YBaaCuaOy 
[ill, La 2 - :c Ba a; Cu04 H3 and HgBa 2 Cu0 4 +5 H3 we 
know that at optimal doping the saddle point is pinned to 
the Fermi level. Since in Hgl201 the vHS reaches Ep at 
approximately 5GPa, we interpret this value as the effec- 
tive "chemical" pressure corresponding to optimal dop- 
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FIG. 4: Densities of states at the Fermi level No in states per 
eV and unit cell in Hgl201, Hgl212, Hgl223, and Hgl234 for 
pressure values of 0, 5, 10 and 15 GPa. 



ing in the sense that both the optimal doping and this 
pressure lead to similar changes in the electronic band- 
structure in the vicinity of the Fermi level, which are im- 
portant for superconductivity. The changes of T c due to 
pressure occur on top of the doping-induced T c , i.e. are of 
the order of 10-20%. The same is observed in Nq when in- 
creasing pressure above 5 GPa as shown in Fig. From 
a comparison of the pressure-induced change of T c with 
the change of No we can estimate A = ANq/AT c x T c /Nq 
being of the order of one. Taking into account a possible 
pressure enhancement of ujb can only increase this esti- 
mated value. Therefore, we can rule out the BCS-like 
weak coupling approaches. 

Second, we want to comment on the van Hove scenario. 
We find that it is not the driving force for boosting T c in 
the sense that it never shows up as a sharp peak in the 
density of states at the Fermi level. At ambient pressure, 
there exists only a saddle point in the band structure, i.e. 
a rather broad feature in the density of states which is 
further considerably broadened when pressure is applied. 
Thus, the effect of a non-uniform density of states near 
the Fermi level, which was supposed to be important to 
enhance T c |l4j in comparison with standard BCS theory, 
is not applicable in this case. Nevertheless, this saddle 
point plays an important role as it is pinned to the Fermi 
level at optimal doping. The pressure effect on T c can 
be traced back to the above discussed pressure-induced z 
dispersion and the corresponding changes in the hopping 
matrix elements. 

Summarizing our results, we find the hole content in 
the Cu02 planes to increase as a function of pressure by 
approximately 0.05 electrons going from ambient pres- 
sure to 15 GPa. These changes in the hole concentration 



are accompanied by a considerable gain in the density 
of states at the Fermi level arising due to the movement 
of the saddle point up to and even through the Fermi 
level. The latter effect is caused by an increase in the 
z axis dispersion due to pressure. We find the effect of 
optimal doping to be equivalent to applying a pressure 
of approximately 5GPa. On top of this, the increase in 
the density of states is of the order of 10% in parallel 
to the enhancement of T c which is of the same order of 
magnitude. From this comparison we conclude that su- 
perconductivity in this cuprate family is described by the 
intermediate to moderately strong coupling regime, i.e. 
A being of the order of 1. The saddle point, often referred 
to as a vHS, does not show up as a peak in the density 
of states making models based on a non-uniform density 
of states being not applicable to high T c cuprates. 



Acknowledgments 

We appreciate support from the Austrian Science Fund, 
projects P13430, P14004, and M591, and by NSF Grant 
No. DMR-02-05125. Most of the calculations were done 
at the Materials Simulation Center, a Penn-State MR- 
SEC and MRI facility. C.A.D thanks E. Antipov for 
valuable discussions. 



[1] 

[2] 

[3] 

[4 
[5 

[6 

[7: 

[8 
[9 
[10 

[11 
[12 
[13 

[14 



S.N. Putilin, E.V. Antipov, O. Chmaissem, and M. 
Marezio, Nature 362, 226 (1993). 

L. Gao, Y.Y. Xue, F. Chen, Q. Xiong, R.L. Meng, D. 
Ramirez, and C. W. Chu, Phys. Rev. B 50, 4260 (1994). 
C.W. Chu, L. Gao, F. Chen, Z.J. Huang, R.L. Meng, 
and Y.Y. Xue, Nature 365, 323 (1993). 
J. Labbe and J. Bok, Europhys. Lett. 3, 1225 (1987). 
P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, 
and J. Luitz, WIEN2k, An Augmented Plane Wave + 
Local Orbitals Program for Calculating Crystal Properties 
(Technische Universitat, Vienna, 2001). 
David J. Singh, Planewaves, Pseudopotentials and 
the LAPW Method, Kluwer Academic Publishers, 
Boston/Dordrecht/London, 1994. 

C. Ambrosch-Draxl, Physica Scripta T 109, 48 (2004). 
T. Thonhauser, H. Auer, E. Ya. Sherman, and C. 
Ambrosch-Draxl, Phys. Rev. B 69, 104508 (2004). 
E.V. Antipov, A.M. Abakumov, and S.N. Putilin, Su- 
percond. Sci. Technol. 15, 31 (2002). 

C. Ambrosch-Draxl, P. Siile, H. Auer, and E. Ya. Sher- 
man, Phys. Rev. B 67, 100505 (2003). 

R. Kouba, C. Ambrosch-Draxl, and B. Zangger Phys. 
Rev. B 60, 9321 (1999). 

T. Thonhauser and C. Ambrosch-Draxl, Phys. Rev. B 
67, 134508 (2003). 

D. L. Novikov, M.I. Katsnelson, Jaejun Yu, A.V. Post- 
nikov, A.J. Freeman, Phys. Rev. B 54, 1313 (1996). 

M. L. Horbach, Phys. Rev. B 47, 15312 (1993). 



